Spin Currents Induced by Nonuniform Rashba-Type Spin-Orbit Field 
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We study the spin relcixation torque in nonmagnetic or ferromagnetic metals with nonuniform 
spin-orbit coupling within the Keldysh Green's function formalism. In non-magnet, the relaxation 
torque is shown to arise when the spin-orbit coupling is not uniform. In the absence of an external 
field, the spin current induced by the relaxation torque is proportional to the vector chirality of 
Rashba-type spin-orbit field (RSOF). In the presence of an external field, on the other hand, spin 
relaxation torque arises from the coupling of the external field and vector chirality of RSOF. Our 
result indicates that spin-sink or source effects are controlled by designing RSOF in junctions. 
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I. INTRODUCTION 



Generation of spin currents is of essential importance in spintronics. So far, several methods have been known to 
generate spin currents, namely, spin injection from ferromagnetic metals, the spin Hall effect the spin pump- 
ing effect [3, and the spin injection by the spin accumulation (spin chemical potential) [aJM- It was recently 
demonstrated that thermal generation of spin currents is also possible (the spin Seebeck effect) [7[. 

In solids, the spin density (s) and spin current (jig) satisfy the continuity equation with a source or sink term, 7"; 

s + v-i, -r. (1) 

The right-hand side is the spin relaxation torque. In metals, it arises mainly from the spin-orbit interaction. In 
the case of uniform and dynamic magnetization, the relaxation torque reduces to the Gilbert damping torque 0, Q . 
In the context of the spin-transfer torque in inhomogeneous spin structures, the spin relaxation torque leads to a 
non-adiabatic torque {(3 term) ^, . In the case of electric spin injection from ferromagnetic metal to non-magnetic 
metal, the spin relaxation torque was shown to be proportional to V ■ E [E is the applied electric field) and the spin 
injection effect is explained based on the behavior of spin relaxation torque at the interface [lo| . The spin relaxation 
torque thus plays an essential role in spin current generation and dissipation. 

To detect the spin current electrically, the inverse spin Hall effect is widely used by attaching metals with strong 
spin-orbit interaction, such as Pt In the non-local spin injection experiment, Kimura et al. discussed on a 

phenomenological ground that the spin current in Cu is absorbed by Pt due to the strong spin relaxation in Pt (spin 
sink effect), and that the absorbed spin current is converted into electric voltage by the inverse spin Hall effect [l^ . 
What is essential in the spin sink effect is a spatial inhomogeneity of the spin-orbit interaction, namely, the gradient 
of the spin-orbit interaction of the spin relaxation torque. In this paper, we investigate this possibility by microscopic 
calculation. 



II. MODEL 

We consider the Hamiltonian with Rashba spin-orbit interaction cou pled with an electromagnetic field. RSOF 
exists when the inversion symmetry is broken at the interface or surface jl3l . The total Hamiltonian T-L is composed 
of non-perturbative part + + ^imp and perturbative part + ^so ^~ ^cm- Each term is 

^o=E J drcUr,t)(^-^\/^-eF^c,{r,t), (2) 

J2 [ dr aMcl{r,t)c,ir,t), (3) 

f dr u(r)4(r,i)c,(r,i), (4) 

a=±l 

h J2 h^'^ -E^so(r-) •4(r,t)(V X 

a,/3=±l •' 

J2 f dr Esoir)-ci{r,t){A{r,t) xa^fs)cp{r,t), (6) 

Hom = ^^ j dr A{r,t)-c\{r,t)Vc„{r,t), (7) 

where Ca{r^t) represents the annihilation operator of a conduction electron, <Tq,^ represents Pauli matrices, M is the 
magnetization, Eso{r) is RSOF and A{r,t) is a vector potential for the electromagnetic field. u{r) is the random 
impurity potential and averaging over the impurity position is carried out as {u{r)u{r'))i^-p = UQnijnpS{r — r'), where 
Mo and riimp are the strength of the impurity and the impurity concentration respectively. The RSOF is assumed to 
arise from the gradient of a scalar potential and thus to satisfy V x Eso = 0- In the following, we calculate the spin 
torque in both systems, that is, non-magnets and ferromagnets, where we adapt M = for non-magnets and M ^ 
for ferromagnets. 

The spin density in the continuity equation is defined as s" = '^pj{c'fjO''^jC^), where f3 and 7 represent the spin 



Hm — - 

(T = ±l 



a,/3=±l 
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indices, and a spin current polarized in the a direction flowing in the i direction is defined by 

{Jf)^ ^ ^^{cHr,t)a"d,c{r,t)) - ^A{r,t){cHr,t)a'^c{r,t))-J2e,^.El,o{r){cHr,t)c{r,t)). (8) 

j 

The spin relaxation torque is derived by taking quantum and statistical averages of the Heisenberg's equation of 
motion for the spin operator Q. Its a component reads 

= * E E ^vkeck^E'soir) (c^ (r, t)a^d,c{r, t)) + y E E ^^Jk^o.k^El,o{r)A, (r, t) {c\r, t)a^c(r, t)), (9) 
where tijk stands for the Levi-Civita symbol. 

A. Equilibrium spin currents 

By calculating the spin relaxation torque perturbatively within the Keldysh Green's function formulation, we access 
the issue of spin source and sink. The contribution of spin relaxation torque which is i-th order in the RSOF and j-th 
order in the vector potential is denoted as j)- First, we deal with the case without the vector potential, i.e., 7^"o)- 
The first order contribution with respect to the RSOF, T^^y vanishes. The second order contribution is calculated 

as 

"^2,0) ~*2 EE E ^ijk^ak-y^lmn J 
i,j~k 7 l,m.n 

E E E e*^^+^^^'''(2fc - Pi), (2fc + Pi)™i?^o(p)4o(Pi)tr [fT^5fc.-^".9;i-p,.-]< , (10) 

k p pi 

where gk^^j denotes Green's function of free conduction electrons in the wave-number space. Trace in the spin space 
is represented by tr and [ denotes the lesser component. 

1. Non-magnets 

We first consider non-magnets, where g^j^ — gk.u^ (1 is the identity matrix of spin space). By using the rotational 
symmetry, the angular average of the wave vectors is calculated as (kikjkikm) = j^k^{5ijSim + SuSjm + SimSji). We 
see that the leading contribution contains the second order derivative of the RSOF. The torque is calculated as 

^-0) = 4 (S)7 ^ E e'^^^-^-E^' [(aL)' - i9L)'] (Eso X plEsoiP^r , (H) 

where g^^ represents the advanced Green's function and gj^\^ is the retarded one. The result is 

27r?i v 

where v represents the density of states. As we see, the spin relaxation torque of eq. (jl4p is zero when the RSOF is 
spatially uniform. The relaxation torque of eq. (jl4l) generates, when the system is static, the spin current given by 

{jfh^^-iEsoxid.EsoT- (13) 
15 ep 

This implies that vector chirality of the RSOF, Eso x {diEgo), produces a spin current. 

2. Ferromagnets 

In the case of spin-polarized conduction electrons, i.e., in ferromagnets, the spin relaxation torque is affected by 
spin polarization. Nevertheless, the contributions of the zeroth and first order derivative of the RSOF vanish in the 
same way as in unpolarized case, resulting in 

r,%) -^{tF + T^) i^^o X V^EsoT ■ (14) 



■7^2,0) = ^77 (^soxV^£;so)", (12) 
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(a) 




FIG. 1. Diagrammatic representation of the contribution of the spin relaxation torque. Contribution (a) is the second order 
contribution in the RSOF. Contribution (b) and (c) are the first order in the RSOF (dotted hues) and hnear in the external 
field (wavy lines). 



Here, v„ {a — ±) is the density of states for spin-split bands and = ep ~ crM. Under equilibrium condition, 
whether the spin of conduction electrons is polarized or not is not so important for the qualitative behavior of the 
spin relaxation torque. 



B. Non-equilibrium spin currents 



We study in this subsection the contribution linear in the external field. 



1. Ferromagnets 



First, let us discuss the spin-polarized conduction electrons. In this case, the leading contribution of spin relaxation 
is the first order with respect to RSOF as shown in Figs. [1] (b) and[T] (c). The spin relaxation torque then reads 



27r / 27r 



''Y.Y.Y. e'^^+''^"(2fc - P)j (2fc + q)i Mq, ^)El,o{p)ty K.9fc.c.gfc-q,.-n] 

k p q 



(15) 



Using (kikj) — ifc^Jy as a result of rotation symmetry, the leading contribution of the torque is 



7(1,1) - * 



2eh f dQ 
3m / 2tt 



{Eso{p) X A{q,n)) X e. 



k,<7 



(16) 



Here, denotes the unit vector of z direction. We obtain 



{Esq X iJom) X 



(17) 



where Tct((T — ±) is the relaxation time for spin-split bands. As we see, the spin relaxation torque as a linear response 
of an external field arises even in the uniform Rashba field case. This spin relaxation may be the dominant mechanism 
of spin sinks in the junction of a ferromagnet and a non-magnet where Sem is applied along the ferromagnet and the 
interface Rashba exists perpendicular to Ecm.. We note that contributions containing derivatives such as in Eq. (12) 
arise also in ferromagnets as so small corrections to 77" 
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FIG. 2. Diagrammatic representation of the contribution of the second order in the RSOF and hnear in the external field. 



2. Non-magnets 

If we assume that A(r, t) — A{t) for simplicity, the dominant torque in the spin-unpolarized case is diagrammatically 
given by (a) and (b) in Fig. [2] and is described as 

i.j.k 7 1,171,11 o k p Pi,q 

X [{2k -pi - q)j{2k+pi)„ii2k - 2pi + q)o £;so(p)^so(Pi)^o(9' 9k,uj(T'^ gk-pi,i^gk-pi-q.i^-n]'^ 

+{2k -pi- q)j{2k - 2q + pi),n{2k + q)o £:^o(p)^^So(Pi)^o(9> [<^'^gk,u:9k-q,LJ-nO-"'gk-q-pi.LJ-n]^l^8) 

By noting V x Eso = and V x Ecm — 0, the leading contribution of the torque turns out to be the forth order in 
wave number and the first order in the derivative of RSOF, resulting in 



P,Pi \ k 

(Esoip) X {A{n)-p,)EsoiPiT + A{n) ■ (Esoip) x d^EsoiPi))]- (19) 
Here, gj^ = gj^^^^o- We obtain 

TTg 7/ r 

7^2 1) = 77— (^^SO X {E,n. ■ V)Esor + ^^cm • (-Eso X daEso) ■ (20) 

The structure of the first term of the right hand side of eq. (|20|) is the same as that of so-called f3 term in the current- 
induced torque, —13{S x (j • V)S'), if we replace the localized spin S by Eso- The second term represents coupling of 
the electromagnetic field and the vector chirality of RSOF. 

Let us look in detail into the spin relaxation torque of eq. (j20p for a junction of two nonmagnetic metals shown 
in Fig. 3 (a). As seen from eq. (l20l) . the RSOF at the interface needs to form a vector chirality {Ego x dEso) to 
induce finite spin relaxation torque. Such configuration is expected to arise generally when the junction is of finite 
area. Since the Rashba effect is due to the structure-inversion-symmetry breaking of a quantum well and the RSOF 
arises from the gradient of a scalar potential, the RSOF will behave near the edges as shown in Fig. 3 (b). Let us 
consider a square junction area in the xy plane with an applied electric field (or current) along the y direction. Then 
both the first and the second terms of the right-hand side of eq. (f20|) have only the component with a = x. The sign 
of the relaxation torque is negative along the x direction and positive along y direction, resulting in the torque density 
distribution of Figs. 3 (e) and 3 (f) with the spin sink along the x axis and spin source along the y axis. In the static 
case, the spin current generated by the present torque density is polarized in the x direction and satisfies V • — . 
In the analogy with the Gauss's law of the classical electromagnetism, the generated-spin-current distribution is as 
depicted in Figs. 3 (c), 3 (d) and 3 (g). Since the induced torque comes from the vector chirality of the RSOF, the 
spin-current distribution generated by the spin relaxation torque is the same for the opposite sides of the system. 



6 




FIG. 3. (Color online) Schematic figures of spin currents induced by tlie spin relaxation torque of eq. (|20p . RSOF is induced 
in the interface between the upper non-magnet and the lower one. Yellow arrows denote the RSOF (Esq), a red arrow denotes 
an external field (Ecui) and purple arrows denote spin currents (jT)- Orange (blue) region indicates plus (minus) values of 
x-component spin relaxation torque, (a): A square junction in the plane, (b): The spatial distribution of the RSOF near the 
edge, (c) and (d): Configuration along the x axis and y axis respectively, (e) and (f): The spatial distribution of a;-component 
spin relaxation torque, (g): The spin-currents distribution induced by the spin relaxation torque. 



which is different from that induced by the spin Hall effect. Experimental observation of the local spin current would 
be interesting. 

As we have discussed, the spin continuity equation is useful to understand the spin current generation and absorption 
based on the spin relaxation torque distribution. One should note, however, that the result of the spin current 
determined this way has an ambiguity given by a rotation of some vector, V x C (C is any vector field). It would 
be important to investigate other equations for the spin current defining its rotation. (In the case of static classical 
electromagnetism, the electric field is constrained by the rotation free condition, V x E = 0.) 



III. SUMMARY 



In summary, we presented a microscopic theory of spin relaxation torque in nonmagnetic and ferromagnetic metals 
with nonuniform spin-orbit coupling within the Keldysh Green's function formalism. We found that the spin relaxation 
torque depends quantitatively on whether an external electric field is applied or not and on the spin polarization of 
the carrier. In the case of non-magnets, there exist equilibrium spin currents induced by vector chirality of RSOF 
in the absence of an external field. In the presence of an external field, on the other hand, spin relaxation torque 
arises from the coupling of the electric field and vector chirality of RSOF. In ferromagnets, in contrast, uniform RSOF 
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results in the spin relaxation torque. Our result indicates that control of RSOF will make possible the spin-currents 
control without magnetism or magnetic field. 
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